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Introduction 

With increasing use of ceramic materials in high temperature structural applications such as advanced heat engine 
components, the need arises to accurately predict thermomechanical behavior (Fig. 1). This paper, which is the first of 
two in a series, will focus on inelastic deformation behavior associated with these service conditions by providing an 
overview of a viscoplastic constitutive model that accounts for time-dependent hereditary material deformation (e.g., 
creep, stress relaxation, etc.) in monolithic structural ceramics. 

Early work in the field of metal plasticity indicated that inelastic deformations are essentially unaffected by 
hydrostatic stress. This is not the case, however, for ceramic-based material systems, unless the ceramic is fully dense. 
The theory presented here allows for fully dense material behavior as a limiting case. In addition, ceramic materials 
exhibit different time-dependent behavior in tension and compression. Thus, inelastic deformation models for ceramics 
must be constructed in a fashion that admits both sensitivity to hydrostatic stress and differing behavior in tension and 
compression. A number of constitutive theories for materials that exhibit sensitivity to the hydrostatic component of 
stress have been proposed that characterize deformation using time-independent classical plasticity as a foundation. 
However, none of these theories allow different behavior in tension and compression. In addition, these theories are 
somewhat lacking in that they are unable to capture creep, relaxation, and rate-sensitive phenomena exhibited by 
ceramic materials at high temperature. 

When subjected to elevated service temperatures, ceramic materials exhibit complex thermomechanical behavior 
that is inherently time-dependent, and hereditary in the sense that current behavior depends not only on current 
conditions, but also on thermo-mechanical history. The objective of this work (Fig. 2) is to present the formulation of a 
macroscopic continuum theory that captures these time-dependent phenomena. Specifically, the overview contained in 
this paper focuses on the multiaxial derivation of the constitutive model, and examines the scalar threshold function and 
its attending geometrical implications. For complete details of the model, the reader is directed to the recent publication 
by Janosik and Duffy (ref. 1). 

Flow Potential 

Using continuum principles of engineering mechanics, the complete viscoplastic theory (Fig. 3) is derived from a 
scalar dissipative potential function, identified here as fi (Fig. 4), first proposed by Robinson (ref. 2) for a J2 model, 
and later extended to sintered powder metals by Duffy (ref. 3). The specific form adopted for the flow potential is an 

’ NASA Resident Research Associate at Lewis Research Center. 


1 


Paper 15 



integral format proposed by Robinson (ref. 2) that has similar geometrical inteipretations (e.g., convexity and 
normality) as the yield function encountered in classical plasticity. This isothermal formulation includes parameters for 
viscosity, hardening, recovery, and unitless stress exponents. Also included is the octahedral threshold shear stress, K , 
which is generally considered a scalar state variable that accounts for isotropic hardening (or softening). However, since 
isotropic hardening is often negligible at high homologous temperatures (> 0.5), to a first approximation K is taken to 
be a constant for metals. This assumption is adopted in the present work regarding ceramic materials. Specific details 
regarding the experimental test matrix needed to characterize these parameters will appear in a second article. 

Threshold Function 

The specific formulation used here for the threshold function, F ( a component of the flow potential function), was 
originally proposed by Wiliam and Wamke (ref. 4) in order to formulate constitutive equations for time-independent 
classical plasticity behavior observed in cement and unreinforced concrete. Wiliam and Wamke (ref. 4) proposed a 
yield criterion for concrete that admits a dependence on the hydrostatic component of stress and explicitly allows 
different material responses in tension and compression (Fig. 5). Several formulations of their model exist, i.e., a three- 
parameter formulation and a five-parameter formulation. For simplicity, the work presented here builds on the three- 
parameter formulation shown in Fig. 6. 

The Wiliam- Wamke criterion uses stress invariants to define the functional dependence on the Cauchy stress (cry) 
and internal state variable (ay). The invariants Tj and admit a sensitivity to hydrostatic stress, and account for 
different behavior in tension and compression (since this invariant changes sign when the direction of a stress 
component is reversed), respectively. Note that a threshold flow stress is similar in nature to a yield stress in classical 
plasticity. The specific details involved in deriving the final form of the function F can be found in Wiliam and Wamke 
(ref. 4). A similar functional form is adopted for the scalar state function G. This formulation assumes a threshold does 
not exist for the scalar function G, and follows the framework of previously proposed constitutive models based on 
Robinson’s (ref. 2) viscoplastic law. 

For the Willam-Wamke three-parameter formulation, the model parameters include oy, the tensile uniaxial 
threshold stress, <j c , the compressive uniaxial threshold stress, and ofr c , the equal biaxial compressive threshold stress 
(Fig. 7). The Willam-Wamke model yields a flow surface in the shape of a pyramid with a triangular base in the Haigh- 
Westergaard stress space, as depicted in Fig. 5. As a reference, typical J 2 plasticity models have yield surfaces that are 
right circular cylinders in the Haigh-Westergaard stress space. 

Flow Surfaces - Interpretation 

As in Robinson’s original theory, the current model is closely tied to the concepts of a potential function and 
normality. It is this potential-normality structure that provides a consistent framework. According to the stability 
postulate of Drucker (ref. 5), the concepts of normality and convexity are important requirements which must be 
imposed on the development of a flow or yield surface. Constitutive relationships developed on the basis of these 
requirements assure that the inelastic boundary-value problem is well posed, and solutions obtained are unique. 
Experimental work by Robinson and Ellis (ref. 6) has demonstrated the validity of the potential-normality structure 
relative to an isotropic J 2 alloy (i.e., type 316 stainless steel). With this structure, the direction of the inelastic strain rate 
vector for each stress point on a given surface is directed normal to the flow surface F= constant , as illustrated in Fig. 8. 
Without experimental evidence to the contrary, it is postulated that this structure is similarly valid for isotropic 
monolithic ceramic materials. The convexity of the proposed flow surface assures stable material behavior, i.e., 
positive dissipation of inelastic work, which is based on thermodynamic principles. The convexity requirement (Fig. 8) 
also implies that level surfaces of a function are closed surfaces, since an open region of the flow surface allows the 
existence of a load path along which failure will never occur. Finally, the Willam-Wamke flow criterion (and the 
constitutive theory presented herein) degenerates to simpler models (e.g., the two-parameter Drucker-Prager and the 
one-parameter Von Mises) under special limiting conditions (Fig. 9). 

Constitutive Model 

Constitutive equations formulated for the flow law (strain rate) and evolutionary law are given in Fig. 10. These 
relationships employ stress invariants to define the functional dependence on the Cauchy stress and a tensorial state 
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variable. The potential nature of Q is exhibited by the manner in which the flow and evolutionary laws are derived. The 
flow law, sjj, is derived from i?by taking the partial derivative with respect to die applied stress. The adoption of a 
flow potential and the concept of normality, as expressed in the flow law, were introduced by Rice (ref. 7). In his work 
this relationship was established using thermodynamic arguments. The authors wish to point out that this concept holds 
for each individual inelastic state. The evolutionary law, a,y, is similarly derived from the flow potential. The rate of 
change of the internal stress is expressed by taking the partial derivative of /2with respect to the internal stress, and 
multiplying by a scalar hardening function dependent on the inelastic state variable (i.e., the internal stress) only. Using 
arguments similar to Rice's, Ponter and Leckie (ref. 8) have demonstrated the appropriateness of this type of 
evolutionary law. 

Example 

A preliminary example was performed to illustrate some of the capabilities of the multiaxial constitutive model 
developed herein, and to compare the results with those obtained utilizing a J 2 model. In addition, the effects of 
varying the Wiliam- Wamke parameters were demonstrated. Figs. 11-13 depict the applied loading condition, 
approach, and observed results. 

A second article will examine specific time-dependent stress-strain behavior that can be modeled with the constitutive 
relationship presented in this article. No attempt is made here to assess the accuracy of the model in comparison to 
experiment. A quantitative assessment is reserved for a later date, after the material constants have been suitably 
characterized for a specific ceramic material. 

Summarv/Conclusion 

The overview presented in this paper is intended to provide a qualitative assessment of the capabilities of this 
viscoplastic model in capturing the complex thermomechanical behavior exhibited by ceramic materials at elevated 
service temperatures. Constitutive equations for the flow law (strain rate) and evolutionary law have been formulated 
based on a threshold function which exhibits a sensitivity to hydrostatic stress and allows different behavior in tension 
and compression. Further, inelastic deformation is treated as inherently time-dependent. A rate of inelastic strain is 
associated with every state of stress. As a result, creep, stress relaxation, and rate sensitivity are phenomena resulting 
from applied boundary conditions and are not treated separately in an ad hoc fashion. Incorporating this model into a 
non-linear finite element code would provide industry the means to numerically simulate the inherently time-dependent 
and hereditary phenomena exhibited by these materials in service. 

References 

1 . Janosik, L.A.; Duffy, Stephen F.: "A Viscoplastic Constitutive Theory for Monolithic Ceramics - 1,” Transactions 
of the ASME Journal of Engineering for Gas Turbines and Power, Paper No. 96-GT-368. Presented at the ASME 
International Gas Turbine Congress, Exposition, and Users’ Symposium, Birmingham, UK, June 10-13, 1996. 

2. Robinson, D.N.: "A Unified Creep-Plasticity Model for Structural Metals at High Temperature," ORNL/TM 5969, 
1978. 

3. Duffy, S.F.: "A Unified Inelastic Constitutive Theory for Sintered Powder Metals," Mechanics of Materials , Vol. 

7, 1988, pp. 245-254. 

4. Wiliam, K.J.; and Wamke, E.P.: "Constitutive Model for the Triaxial Behaviour of Concrete, Int. Assoc. Bridge 
Struct. Eng. Proc., Vol. 19, 1975, pp. 1-30. 

5. Drucker, D C.: "A Definition of Stable Inelastic Material", Journal of Applied Mechanics, Vol. 26; ASME 
Transactions, Vol. 81, 1959, pp. 101-106. 

6. Robinson, D.N.; and Ellis, J.R.: "Experimental Determination of Flow Potential Surfaces supporting a Multiaxial 
Formulation of Viscoplasticity," Proceedings, 5th International Seminar on Inelastic Analysis and Life Prediction 
in High Temperature Environments, Paris, France, 1985. 

7. Rice, J.R.: "On the Structure of Stress-Strain Relations for Time-Dependent Plastic Deformation in Metals," 
Journal of Applied Mechanics, 37, 728, 1970. 

8. Ponter, A.R.S.; and Leckie, F.A.: "Constitutive Relationships for Time-Dependent Deformation of Metals," ASME 
Journal of Engineering Materials and Technology, Vol. 98, 1976. 


3 


Paper 15 



CERAMIC LIFE PREDICTION ANALYSIS 


POTENTIAL 
FAILURE MODES 
IN HEAT ENGINE 
APPLICATIONS: 


Time-Independent Failure Modes: 

• Fast fracture (tension, compression, shear) 

• Buckling 

Time-Dependent Failure Modes: 

• Static, dynamic, and cyclic fatigue 

• Creep crack growth 

• | Creep deformation response | 

• Stress corrosion and oxidation 

• Impact and contact loading response 


Current Capabilities: 

Needed Capabilities: 

✓Fast-fracture 


■ ■ J - 

•| Creep deformation | 

✓Subcritical Crack Growth 


• Creep rupture 

High loading 


Extended duration loading 

Low to intermediate temperatures 


Elevated temperatures 

(-2000 °F) 


(-2400 °F) 


Fig. I 


OBJECTIVE 

To present a multiaxial continuum theory which 
accounts for time-dependent hereditary material 
deformation behavior of isotropic structural ceramics. 


OUTLINE: 


• Introduction to Viscoplasticity 

• Flow Potential 

• Threshold Function 

• Flow Surfaces-lnterpretation 

• Constitutive Model 

• Example 

• Conclusions/Enhancements 


Fig. 2 
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VISCOPLASTICITY 


Phenomenological 

Scalar 


Viscoplastic 

Approach: 

Potential 

Continuum 

Constitutive 


Function 

principles 

Model 


(Considers effects at the macrostractural level) 



Fig. 3 


FLOW POTENTIAL FUNCTION 


General form: 

Q = a(o ij ,a ij ,T) 

Assume isothermal 

conditions 


Q = Q fcj.ttij) 


applied/ \ 

stress ' x 


internal 

state 

variable 


Specific 

form: 


(Robinson, 1978) 


where fi, R,H,n,m, and K are material constants 
characterizing viscoplasticity 



where: = Sij -ay Sg = a i} - (1/3) <5, y 

Vij = a ij = a ij ~ (1/3)tt « 5 U 
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THRESHOLD FUNCTION 


• Originally proposed to formulate Willam-Wamke Model 

constitutive equations for time- 
independent classical plasticity behavior 
in cement and unreinforced concrete 

• Three-parameter model 

• Assumes material is isotropic 

• Allows different behavior in tension & 
compression through the invariant J 3 

• Admits dependence on the hydrostatic 
stress through the invariant h 




when 1 } = (1/3) 0u 

J 2 - (1/2) Sy Sji 
J 3 = (1/3) Sfifrt 

Fig- 5 



THRESHOLD FUNCTION 


Inelastic deformation occurs for 
stress states where 1 




Willam-Wamke 
3-Parameter Formulation: 


Define: 


'1 = Vu 


Jo = 


'll z » x * 




I,. 

31 y jk ki 


F J 2 . J 3 ) = j- 

u c 


9 f 1 

1 

1 

, i 

3 PO c | 


(Wiliam and Wamke, 1975) 


The dependence on J 3 is introduced through the angle of similitude 9 

(^F)h 

Note that similar relationships exist for scalar function G. 


cos(3e) =tjr. 3/2 


Fig. 6 
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THRESHOLD PARAMETERS 



Fig. 7 

FLOW SURFACES 

Normality Condition: 

The direction of the inelastic strain rate vector for each stress point on a given surface is 
directed normal to the flow surface F=constant. 


• Flow surfaces eventually cluster forming a limiting surface. 

• Implies large changes in elastic strain rate for only small stress changes. 

Fig. 8a 
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FLOW SURFACES 

Convexity Condition: 


Willam-Warnke: OJ>r c <r t <r c 



• Assures stable material behavior (i.e., positive dissipation of inelastic work). 

• Implies that level surfaces of a function are closed surfaces. 

Fig. 8b 


FLOW SURFACES 

(Degeneration to “simpler” models) 

Limiting Conditions: 



Wiliam- Warnke 
(‘Triangle” in 7t-plane) 
3 Parameters 
r t, r c , p 


Drucker-Prager 
(Circle in n-plane) 
2 Parameters 

T o>P 

l 

Von Mises 
(Circle in Ji-plane) 

1 Parameter 

r 0 

r c -r t = r 0 

Fig. 9 

r c = r t = r 0 

P —>oo 
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CONSTITUTIVE MODEL 

The potential nature of £2 is exhibited in the extended normality structure 


Flow Law: 


i M/rr xxi rr * 

d Q 

Eij " dCij 


&ij “ ^0 

C\Sij + C 2 I,y+ C 3 

L _ 27 2 %11 

3 n 


[Rice (1967), (1971)] 


Limiting Conditions: ^ ^ J 2 Model 

T c ~ r t~r 0 and p Ey = Co[C 2 £y] 


Evolutionary ha w: 


d Q 

“y “ _h day 


dij = h- 

£lj C 4 

C]5y + C 5 a ij + C 6 

^ j ^ 

3 J 


[Ponter & Leckie (1976) > 



Ponter (1976)] where h is a hardening function dependent on the internal stress 


Limiting Conditions: 
r c = r ( = r 0 andp-> c 


a 


Fig. 10 


y 


J 2 Model 

= h {£y-Q[c 5 «y]} 


Example: Willam-Warnke Model vs. J 2 Model 




10 0 0 
000 
000 


ksi 


Willam-Warnke Model: 


Subcase [X\ 

Subcase B 


Subcase C 

o, - 0.2 

0 , = 0.2 

a, = 0.5 

a c = 2.0 

c £ = 2.0 

o c = 1.0 

o be =2.32 

it 

*0 

II 

0 " 


J 2 Model: 


o, = 1.732 
c c = 1.732 
G bc = 1.732 


Fig. II 


APPROACH 


Viscoplasticity 

parameters 

Pf W; n, tn, b 

Willam-Warnke 

parameters 

(<%, g c , <*bc) or (p,r t ,r c ) 


i 

( Computer Algorithm 
(Constitutive Model) 


• Creep Curve (Strain vs. Time) 

• State Variable vs. Time 

• Flow Surface 
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Willam-Wamke Model vs. J 2 Model 

CUTTING PLANE VIEW 


J2 Model 

• Threshold Parameters 

Ot — Oc— Obc 


Tensile Meridian. 6=0 
Oc 

hhr 

r"~ \ 

/4 

/' _(./ ■ " — ■ ■■ y 

^ 

p — + <*> 

Compressive Meridian, 6=0 



p — * °° 


• Meridians are linear 

• Meridians are parallel 


Willam-Wamke Model 
•Threshold Parameters 


Of ^ (Jc ^ Oftc 

(For typical ceramics, <j,<G c <o bc ) 


Tensile Meridian, 6=0 ^ 


Compressive Meridian, 0= 0 

p •* — 


p = a ‘ a * 

°c «Jbc-<7 t ) 


• Meridians are linear 

• Meridians intersect hydrostatic axis 


Willam-Wamke Model vs. J 2 Model 

J 2 Model 


• Poisson effect (v~0.5) 

• £ll — -2£22 = "2£33 



Fig. 13a 
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Willam-Wamke Model vs. J 2 Model 

Willam-Wamke Model 



Fig. 13b 


SUMMARY 

• Fundamental concepts of viscoplasticity 

• Specific form for the potential function 

• Multiaxial derivation of viscoplastic constitutive model 

• Preliminary applications of constitutive model 


FUTVRE ENDEAVORS 

• Incorporate model into computer algorithm 

• Develop experimental test matrix 

• Identify parameters applicable to ceramics 

• Model phenomena such as creep , stress relaxation, 
strain-rate sensitivity, etc. 

• Develop user-defined subroutines for commercial 
FEA software packages 

Fig. 14 
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